MICROCOPY  RESOLUTION  TEST  CHARI 

NATIONAL  BURiAU  Of  StANdABOs  .-Kbi  A 


one  FILE  COPY  AD-A141  010 


FTD-ID(RS)T-0250-84 


FOREIGN  TECHNOLOGY  DIVISION 


A  MIXED  EXPLICIT-IMPLICIT  ANTIDIFPUSIVE  METHOD  OF 
NAVIER-STOKES  EQUATIONS  FOR  SUPERSONIC  AND 
HYPERSONIC  SEPARATED  FLOWS 

by 

Chang  Han-hsin,  Lu  Lin-sheng,  et  al 


Approved  for  public  release; 
distribution  unlimited. 


FTD  -ID(RS)T-0250-84 


EDITED  TRANSLATION 


FTD-ID(RS)T-0250-84 

MICROFICHE  NR:  FTD -8 4 -C -0004 05 


19  April  1984 


A  MIXED.  EXPLICIT-IMPLICIT  ANTIDIFFUSIVE  METHOD  OF 
NAVIER-STOKES  EQUATIONS  FOR  SUPERSONIC  AND  HYPERSONIC 
SEPARATED  FLOWS 

By:  Chang  Han-hsin,  Lu  Lin-sheng,  et  al 
English  pages:  22 

Source:  Yingyong  Shuxue  he  Lixue,  Vol.  4,  Nr.  1,  1983,  pp. 
53-67 

Country  of  origin:  China 
Translated  by:  SCITRAN 

F33657-81-D-0263 
Requester:  FTD/TQTA 

Approved  for  public  release;  distribution  unlimited. 


FTD  -ID(RS)T-0250-84 


Date '9  April  ^9  84 


GRAPHICS  DISCLAIMER 


All  figures,  graphics,  tables,  equations,  etc.  merged  into  this 
translation  were  extracted  from  the  best  quality  copy  available. 


O  V*.  ■s/  ■.'■  V 


A  MIXED  EXPLICIT-IMPLICIT  ANTIDIFFUSIVE  METHOD  OF  NAVIER-STOKES 
EQUATIONS  FOR  SUPERSONIC  AND  HYPERSONIC  SEPARATED  FLOWS 


/chang  Han-hsin,  ytu  Lin-sheng,  ,Yu  Tse-ch'u,  /tia  Chan-kuei 


ABSTRACT 


In  tiiis  paper,  a  mixed  explicit- implicit  ecliemc  based  oa  aa  aatidiffusire  method  is  used 
to  solre  the  Narier-Stokes  equations  for  the  supersonic  and  hypersonic  separated  flows. 
The  computations  are  performed  for  laminar  and  turbulent  flows  over  the  two-  and  three- 
dimensional  compression  corners.  The  obtained  results  are  compared  with  the  results  of  nu¬ 
merically  computing  NS  equations". >"*'and  these  results  of  the  experiments"".  The 
computations  show  that  the  numerical  scheme  in  this  paper  is  satisfactorr. 


I .  Introduction 

When  a  finite  difference  method  is  used  to  solve  the  problem  of 
separated  flows  of  a  viscous  gas,  the  viscosity  of  the  difference 
scheme  or  the  artificial  viscosity  introduced  must  be  much  smaller 
than  the  actual  physical  viscosity  term  in  the  region  where  the 
effect  of  viscosity  is  important.  Therefore,  it  is  preferable  to 
use  a  scheme  with  a  high  degree  of  accuracy.  If  there  are  shock 
waves  in  the  flow  field,  the  difference  scheme  should  be  able  to 
automatically  capture  the  shock  waves.  Furthermore,  the  scheme  should 
have  a  relatively  long  st2d3le  time  step  length  so  that  computation 
time  will  not  be  very  long  when  computations  are  performed  for  high 
Reynolds  nvimber  laminar  or  turbulent  flows,  despite  the  small  size 
of  the  mesh.  The  method  based  on  a  mixed  explicit- implicit  scheme, 
given  in  [1],  meets  the  first  two  of  the  above  requirements,  and 
has  already  been  used  to  obtain  very  good  results.  However,  computa¬ 
tion  tiine  required  is  rather  long  because  of  the  small  stable  time 
step  lengths.  It  is  the  purpose  of  this  paper  to  seek  a  difference 
scheme  that  can  simultaneously  meet  the  above  three  requirements. 

The  scheme  studied  studied  in  this  paper  is  a  mixed  explicit- 
implicit  scheme.  The  time  split-finite  volume  scheme  is  employed. 
Implicit  expressions  are  used  for  the  difference  operators  along  the 


normal  to  the  wall  in  the  neighborhood  of  the  wall,  while  explicit 
expressions  are  used  for  the  difference  operators  in  other  regions 
and  those  along  other  directions.  One  can  thus  obtain  very  large 
stable  time  step  lengths . 

In  the  implicit  difference  computations,  we  have  used  a  two- 
step  antidiffusive  method  to  increase  the  stability  of  the  computa¬ 
tion  and  the  accuracy  of  the  difference  scheme.  The  first  step  is 
to  adopt  a  first-order  implicit  scheme  containing  a  positive  dissi¬ 
pative  term.  In  the  second  step,  local  explicit  antidiffusion  is 
carried  out,  i.e.,  the  positive  dissipative  term  is  subtracted 
explicitly.  After  these  two  steps,  the  scheme  retains  an  accuracy 
of  the  second  order.  This  type  of  two-step  antidiffusive  method  has 
also  been  applied  to  the  explicit  scheme.  In  ordel:  to  be  able  to 
automatically  capture  the  shock  waves  in  the  computations,  we  have, 
on  the  basis  of  the  first-  and  second-order  schemes  given  by  the 
two-step  method,  established  an  automatically  regulated  mixed  scheme 
Computations  have  shown  that  this  method  is  very  satisfactory. 

This  paper  is  divided  into  four  parts ,• including  the  introduc¬ 
tion.  The  second  part  deals  with  the  two-step  antidiffusive  differ¬ 
ence  method  applied  to  a  model  equation.  In  the  third  part,  the 
explicit- implicit  difference  method  for  solving  the  Navier-Stokes 
(abbreviated  below  as  NS)  equations  is  given.  The  results  of  the 
computation  is  given  in  the  fourth  part  which  also  contains  a  com¬ 
parison  and  analysis  of  the  results. 


II .  Model  equation 


For  the  sake  of  simplicity,  we  consider  the  following  model 
equation: 


^  ,  du 


(2.1) 


In  the  above  equation,  a  is  a  constant.  The  two-step  explicit 
scheme  is  first  established,  then  the  two-step  implicit  scheme  of 
the  antidiffusive  method. 


1.  Two-Step  explicit  scheme  of  the  antidif fusive  method. 


In  the  first  step,  we  adopt  the  windward  scheme,  i.e., 


a>0 

»,*) 

a<0 

(2.2a) 


or  in  a  unified  form. 
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(2. 2b) 
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It  is  not  hard  to  recognize  this  as  the  PycamoB  scheme  for  the 
case  where  u  s  1: 
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It  is  well  known  that  equation  (2.2)  is  a  first-order  scheme. 
Its  viscosity  term  is 

+«<-»)  (2.3) 

where 

(2.3a) 


o.-^( 


The  second  step  consists  of  an  antidiffusion  of  the  first  step, 
i.e.,  the  viscosity  term  (2.3)  is  subtracted  from  u^  .  Hence, 

ur'-irr^— l-Q.Ca-^.-Ear+tt*..)  .  (2.4a) 

Obviously,  equation  (2.4a}  is  accurate  to  the  second  order.  Making 
use  of  equation  (2.2),  we  can  also  write  equation  (2.4a)  as 


aj+tt**'  - 

o>0 

-■^0^-  .TIi)) 

a<0 

(2.4b) 


Combining  equations  (2.2a)  and  (2.4b),  we  can  readily  see  that  this 
is  MacCormack ' s  two-step  explicit  scheme.  This  shows  that  one  can 
arrive  at  the  second-order  MacCormack  scheme  by  starting  out  with 
the  first-order  windward  scheme  and  performing  antidiffusion.  Simi¬ 
larly,  setting  out  from  other  first-order  schemes  containing  positive 
dissipative  terms,  one  can  also  establish  the  corresponding  second- 
order  schemes  by  means  of  antidiffusion. 


2.  Two-step  implicit  scheme  of  the  antidif fusive  method. 


In  implicit  computations,  the  Crank-Nicoson  scheme  is  often 


used: 


-  -■  (2.5) 

This  is  a  second-order  scheme  with  a  growth  factor  of  1..  To  increase 
its  stability  and  to  corroborate  the  advantage  of  the  diagonal  ele¬ 
ments  of  the  matrix  in  the  matrix  chasing,  as  a  first  step,  we  add 
a  positive  implicit  dissipative  term  to  the  right  hand  side  of 
equation  (2,5) .  Thus,  we  obtain  a  first-order  implicit  scheme: 


(2.6) 


Here,  is  the  coefficient  of  dissipation,  u?  can  be  found  from 

equation  (2.6).  Then,  in  the  second  step,  the  dissipative  term 

is  subtracted  partially  explicitly  from  and 

*  .  3 

we  obtain  the  second-order  scheme  below: 

^(••*.-2*? +•?.,)  .  (2.7a) 

Or,  insert  equation  (2.6)  in  equation  (2.7a)  to  obtain 

Or _  _  _  1  ^2. 7b) 

In  order  to.  be  able  to  automatically  capture  shock  waves  in  the  com¬ 
putational  process,  we  have  constructed  on  the  basis  of  equations 
(2.2)  and  (2.4)  or  equations  (2.6)  and  (2.7),  the  automatically 
regulated  mixed  scheme 


In  the  above  equation, 

Q^•9Q^  . 


(2.8) 


(2.9) 


9  is  the  switching  function  in  the  automatic  regulation.  Its  value 
in  the  vicinity  of  the  shock  waves  is  1.  In  regions  where  viscosity 
plays  an  important  role,  its  value  is  zero.  This  shows  that  the 
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scheme  has  accuracy  of  the  first  order  in  the  vicinity  of  shock 
waves,  while  it  has  second-order  accuracy  in  regions  removed  from 
shock  waves-  For  the  expression  for  9,  see  [1].  Computations  have 
shown  that  the  value  of  is  best  taken  to  be  between  0.1  and  0.3. 

In  summary,  we  have  obtained  the  automatically  regulated  anti- 
diffusive  scheme  (2.8)  in  which  can  be  computed  either  explic¬ 

itly  or  implicitly.  It  can  be  shown  by  means  of  a  Fourier  analysis 
that  when  the  explicit  method  is  used,  the  stability  condition  for 
both  equations  (2.2)  and  (2.4)  is  <1.  ,  while  if  the  implicit 

method  is  used,  equations  (2.6)  and  (2.7)  are  unconditionally  stable 


If  the  filtering  function 

•I-*; 


(2.10) 


is  introduced,  then  accurate  to  the  second  order,  equation  (2.6) 
can  be  written  as 


(2.11) 


+yC^?-2*'r'+«'5?) 

Equation  (2.2)  can  be  written  as 

A.a';  (2.12) 

where 

a>0  1 

^  _ _  _  }  (2.13) 

a<0  ‘ 

Equation  (2.8)  becomes 

.  (2.14) 

After  the  above  operation,  equation  (2.14)  becomes  a  second-order 
scheme . 


The  above  study  of  the  model  equation  can  be  extended  to  the 
equation 


Here,  f  is  a  function  of  U.  Proceeding  as  in  equation  (2.6),  we  can 
write  the  first  step  of  the  antidiffusive  implicit  scheme  as 


O’*' 
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after  inserting  it  in  the  above  equation,  we  obtain 


^Au,AU,.,  ( uj:i  _2C/r‘  +C!jTn. 


In  the  above  equation. 


*'  Kdu), 


(2.16) 


(2.17) 


Proceeding  as  in  equation  (2.7a.),  we  obtain  the  second  step  as 


c/r‘-i7r‘  -y(c/u.-2£;’+£/j.,)  . 


(2.18) 


On  the  basis  of  equations  (2.16)  and  (2.18),  one  can  establish  a 
mixed  scheme  as  well  as  introduce  a  filtering  function  and  perform 
filteriiig  computation. 

III.  Explicit-implicit  method  for  solvinq  the  NS  equations 


1.  Basic  equations,  boundary  conditions  and  initial  conditions 

Let  us  now  study  the  diffracted  flows  over  the  two-  and  three- 
dimensional  compression  corners  shown  in  Figures  1  and  2.  As  the 
three-dimensional  corner  possesses  a  plane  of  symmetry,  only  the 
right  half  has  been  shown  in  Figure  1  (as  viewed  along  the  direction 
of  the  oncoming  stream).  Assxime  that  x',  y'  and  z'  are  the  coordi¬ 
nates  of  the  coordinate  systems  given  in  Figures  1  and  2;  t*  denotes 
time;  u'  v'  and  w'  are  the  x'-,  y'-  and  z '-components ,  respectively, 
of  the  velocity  of  the  gas;  p',  p',  T'  and  p'  are,  respectively,  the 
pressure,  density,  temperature  and  coefficient  of  viscosity  of  the 
gas;  y  is  the  adiabatic  index;  +«»''*)•  ;  l'  is  the 

characteristic  length  shown  in  Figures  1  and  2 .  The  following  dimen 
sionless  mesh  quantities  are  introduced: 


iVi 


V  o  .  X/ 


«  =  tn=tu'  lu'  e  —  t!  /»  . 

r— m-mvm'.,  > 

•*  • 

X^X'IL\  y^y'/L\  z^z'/L' 

i^fvJ  ,IU  ' 

where  R  is  the  gas  constant,  and  the  subscript  denotes  the 
oncoming  stream. 


For  the  three-dimensional  diffracted  flow,  the  body  surface 
coordinate  system  C/  n,  ?  is  used.  It  is  related  to  the  Cartesian 
coordinate  system  via  the  following  relations: 

’-V  [  (3.2) 

»)—/(».  y.  *) 

In  the  2d)ove,  z  =  fQ(x,y)  is  the  equation  of  the  body  surface  given 

by  </.)  u:) 

,0,  ■*  t8<P 


jxtgffl.  X>0^>0  ^2  3) 

xtgfij+ytg^fr,  *>05.— *  tg^<V<0  » 

(^,1)— or;  <2)— and 

Here,  tg<t>  =  tgu/tg  .  For  the  significance  of  co  and  <p,  refer  to 
Figure  1.  Therefore,  in  the  n,  C  coordinate  system,  the  NS  equa 
tion  for  the  three-dimensional  laminar  flow  if  a  complete  gas  with 

constant  specific  heat  is 

dU  .  dF  .  dG  ,  9H 

(3.4) 


where 


3.5 


Figure  1.  External  features 
of  the  three-dimensional 
corner  and  its  solution  space 
and  mesh. 

1 —  main  wing  surface; 

2 —  auxiliary  wing  surface; 

3 —  bottom  surface 


L. 


I  u— neca  .  la— 0  x<0 
A  K-i-xtfo  to— o  '>• 


Figure  2 .  External  features 
of  the  two-dimensional 
corner  and  its  solution  space 
and  mesh. 
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For  the  simplified  set  of  NS  equations  [1] 
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H,* 
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For  the  complete  NS  equation. 
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dC 

dw 


dv 


d»7 


(3.9) 


In  equations  (3.6)  and  (3.7),  e  =  0,  1  denote,  respectively,  the 
simplified  and  the  complete  NS  equation.  In  equation  (3.9), 
denotes  the  free  stream  Mach  number  and,  denotes  the 

free  stream  Reynolds  number.  is  Prandtl  number  (0.72  for  air) . 
In  the  expression  for  u,  =  114®K.  In  equation  (3.8a), 

S^w+a/x-f-r/,. 

For  the  two-dimensional  flow  shown  in  Figure  2,  as  the  problem 
is  independent  of  y,  ,/«*•/£ v^p,i»p,§=^6i»' 

Therefore,  the  simplified  and  complete  NS  equations  for  the  two- 
dimensional  case  can  be  obtained  from  equations  (3.4)-(3.9)  . 


Suppose  that  in  the  above  set  of  equations,  u  is  replaced  with 
,  where  u_ ,  the  coefficient  of  viscosity  of  laminar  flow. 


,  where  u_ , 


can  be  cojnputed  from  the  second  to  the  last  equatxon  in  (3.9)  ,  and 
the  coefficient  of  viscosity  of  turbulent  flow,  can  be  computed 
by  means  of  a  zero-equation,  one-equation  or  two-equation  model  (see 
[2]) .  Then,  the  above  set  of  equations  can  also  be  used  to  compute 
turbulent  flow  problems. 

The  regions  shown  in  Figures  1  and  2  are  the  solution  spaces . 

The  surface  of  the  body  is  located  at  the  lower  boundary  of  this 
space.  The  condition  of  absence  of  gliding  flow  is  satisfied  on  the 
wall  surface,  the  temperature  of  which  is  given.  For  the  three- 
dimensional  flow  in  Figure  1,  as  the  corner  affects  the  flow  upstream 
only  in  a  very  limited  region  upstream,  further  upstream  the  flow  is 
two-dimensional.  Therefore,  at  the  cross-section  of  the  entrance, 
with  the  flow  parameters,  given,  the  results  for  a  two-dimensional 
slab  may  be  applied.  At  the  boundary  surface  of  the  exit,  we  consi¬ 
der  two  regions:  on  the  main  wing  surface,  the  flow  does  not  change 
very  much  along  the  C-direction.  Hence,  for  any  flow  parameter  <\>, 
we  have  .  On  the  auxiliary  wing  surface  and  the  bottom  surface, 

the  flow  possesses  conic  characteristics  and  the  physical  quantities 
of  the  flow  are  equal  along  the  radial  lines  originating  from  point 
o.  On  the  plane  of  symmetry,  the  symmetry  requirements  are  satisfied. 
On  the  side  boundaries,  the  flow  does  not  change  much  along  the  n- 
direction,  and  we  have-|^=aO..  On  the  upper  boundary,  because 

it  is  far  from  the  wall  surface.  For  the  two-dimensional  flow  in 
Figure  2,  the  entrance  boundary  may  be  upstream  or  downstreeim  the 
front  edge.  The  flow  parameters  are  given  [1].  At  the  exit  boundary, 

'dtp 

as  it  is  farther  removed  from  the  corner,  we  take  •  On  the 

upper  boundary,  either  single  wave  motion  condition  may 

be  used. 

The  initial  conditions  may  be  taken  to  be  those  for  the  constant 
steady  flow  of  a  nonviscous  gas  around  the  corner,  or  those  for  a 
uniform  flow  field. 

2 .  Mesh  spacing 
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The  mesh  spacing  is  as  shown  in  Figures  1  and  2.  Equally 
spaced  divisions  are  taken  along  the  direction  of  the  stream. 

Smaller  mesh  spacings  are  taken  along  the  ^-direction  near  the  wall, 
while  larger  mesh  spacings  are  taken  in  regions  removed  from  the 
wall.  For  the  three-dimensional  case,  to  better,  describe  the  three- 
dimensional  effect,  smaller  mesh  spacing  is  taken  along  the  horizon¬ 
tal  direction  on  the  auxiliaary  wing  surface. 

3,  Difference  scheme 


According  to  the  time-split  finite  volume  theory,  solving  equa¬ 
tion  (3.4)  is  equivalent  to  solving  the  following  three  equations: 


dU  ,  dF  ^ 

dLi  dG 
dt  *  dn 

^.dH  ^ 
dt  ac  • 


(3.10) 

(3.11) 


(3.12) 


In  regions  near  the  wall,  equation  (3.12)  can  be  further  split  into 

dU  dfik  - 

-jT+ir  “  (3.12a) 


dt  ^  dC  ” 


(3.12b) 


When  solving  equations  (3 ,10 ) -(3 .12) ,  the  mixed  explicit-implicit 
difference  method  is  employed..  The  difference  equations  are 

Here,  ,  is  the  filtering  function,  and  is  related  to  .  .  via 

1,J,K  1,],K 


uT,,,k  -  c/r,,.  +y  Qn(  -2  c/r,*,,.  +  £/?,*>-,.» ) 

-2c;r,/,* -hc/r-,,,..  > 


(3.14a) 


(3.14b) 


(3.14c) 


Refer  to  [1]  for  the  significance  of  Q^,  and  Q^.  The  signi¬ 
ficance  of  Xit  JOii  «C«»‘  are  given  below. 

(1)  are  the  difference  operators  corresponding  to 

equations  (3.10)  and  (3.11),  and  are  defined  at  all  inner  points  of 
the  entire  region  of  computation.  is  the  difference  operator 

corresponding  to  equation  (3.12),  and  is  defined  at  the  inner  points 
of  the  looser  mesh  region  on  the  outside.  These  three  difference 
operators  have  been  obtained  from  a  generalization  of  the  explicit 
scheme  for  the  model  equation  given  in  the  previous  section. 

(2)  jC<p  are  the  difference  operators  of  equations  (3.12a) 
and  (3.12b)  and  are  defined  at  the  inner  points  of  the  tighter  mesh 
region  close  to  the  wall.  The  implicit  scheme  is  used  here.  These 
operators  have  been  obtained  from  a  generalization  of  the  implicit 
antidif fusive  scheme  given  in  the  previous  section.  The  difference 
equation  for  can  be  computed  by  means  of  the  method  of  main  dia¬ 
gonal  matrix  chasing.  As  for  ,  as  equation  (3.12b)  is  parabolic, 
we  may  let  »  0,  and  the  difference  equation  can  be  solved  using 
the  method  of  labeled  quantity  chasing. 

4.  Stability  conditions 

The  stable  time  step  lengths  for 
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In  the  above  equations 
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In  addition,  ^  is  an  inner  point  of  the  solution  space,  and  is 
an  inner  point  of  the  external  looser  mesh  region. 

As  X<»»  are  stable  unconditionally,  the  stable  time  step  ^ 

length  for  the  entire  difference  operation  is 

Af— oin(2A<{»  2Attr  ^0 ,  (3.16) 

Obviously,  this  is  much  longer  them  the  stable  time  step  length  of 
the  explicit  scheme  [1],  [3].  in  particular,  the  former  may  be  made 
1  to  2  orders  of  magnitude  greater  than  the  latter  in  computations 
performed  for  the  turbulent  flow. 

If  we  let  =  0  and  An  -►  <*»  in  equations  (3.15)  and  (3.16)  ,  then 
the  conditions  for  stability  for  the  two-dimensional  flow  can  be 


IV.  Results  of  computation 


Five  examples  have  been  computed  in  this  paper.  The  first 
three  involve  the  laminar  separated  flow  over  the  two-dimensional 
compression  corner.  The  foiirth  is  the  turbulent  separated  flow  of 
the  two-dimensional  compression  corner.  The  fifth  is  the  laminar 
separated  flow  of  the  three-dimensional  compression  corner.  The 
given  conditions  and  the  results  obtained  are  given  in  Figures  3,  4, 
5,  6  and  7.  For  comparison,  we  have  included  in  Figures  3,  4,  5  and 
6  the  results  of  numerical  solutions  of  the  NS  equations  given  in 
[1] ,  [3],  [4]  and  [5],  as  well  as  some  experimental  results  [6],  [7] 
It  can  be  seen  that  the  method  presented  in  this  paper  is  satisfac¬ 
tory.  Compared  with  the  explicit  computational  method,  the  con^juta- 
tion  time  of  our  method  is,  respectively,  4,  13,  9,  20  and  5  times 
shorter  for  these  examples.  Thus,  the  present  method  has  the  advan¬ 
tage  of  saving  machine  time. 


Figure  3(a).  Distribution  of  surface  pressure  in  leuninar  flow 
over  two-dimensional  compression  corner,  with  =  3. 

1 —  solution  of  the  complete  NS  equations  given  in  this  paper; 

2—  solution  of  the  simplified  NS  equations  given  in  this  paper; 

3 —  explicit  solution  of  the  NS  equations  given  in  [1]; 

4 —  explicit  solution  of  the  NS  equations  given  in  [3]. 
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Figure  4  (b)  .  Distribution  of  surface  drag  in  laminar  flow  over 
two-dimensional  compression  corner,  with  =  6.06 

1 —  solution  of  the  NS  equations  given  in  this  paper; 

2—  explicit  solution  of  the  NS  equations  given  in  [4]; 

3 —  adiabatic  wall 


Figure  5(a).  Distribution  of  surface  pressure  in  laminar  flow  over 
two-dimensional  compression  corner,  with  =  14.1. 

1 —  solution  of  the  NS  equations  given  in  this  paper; 

2 —  explicit  solution  of  the  NS  equations  given  in  [3]; 

3 —  explicit  solution  of  the  NS  equations  given  in  [1]; 

4 —  experiiitental  result  [6] 
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Figure  6 (a) .  Distribution  of  surface  pressure  in  turbulent  flow 
over  two-dimensional  compression  corner,  with  =  2.96. 

l--experimental  result  [7];  2 — Shang's  solution  of  the  NS  equations 
(balanced  model)  given  in  [5];  3 — solution  of  the  NS  equations  that 
have  been  simplified  using  a  balanced  algebraic  model  in  this  paper 
4 — solution  of  the  NS  equations  that  have  been  simplified  using  a 
one-equation  model  in  this  paper;  5 — adiabatic  wall 
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Figure  6 (b) .  Distribution  of  surface  pressure  in  turbulent  flow 
over  two-dimensional  compression  corner,  with  =>  2.56. 

1 —  solution  of  the  NS  equations  that  have  been  simplified  using  a 
relaxed  model  in  this  paper; 

2 —  experimental  result  [7] ; 

3 —  adiaUSatic  wall 
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Figure  6 (c) .  Distribution  of  surface  drag  in  turbulent  flow  over 
two-dimensional  coir^jression  corner,  with  »  2.96. 

1 —  solution  of  the  NS  equations  that  have  been  siir5)lified  using  a 
relaxed  algebraic  model  in  this  paper; 

2— -Shang's  solution  of  the  NS  equations  (replaced  mode]  given  in  [5]; 

3—  experimental  result  [7]; 

4—  adiabatic  wall 
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Figxire  6  (d)  .  Distribution  of  surface  drag  in  turbulent  flow  over 
two-dimensional  compression  corner,  with  *  2.96. 

1 —  solution  of  the  NS  equations  that  have  been  simplified  using  a 
balanced  algebraic  model  in  this  paper; 

2 —  solution  of  the  NS  equations  that  have  been  simplified  using  a 
one-equation  model  in  this  paper; 

3—  Shang's  solution  of  the  NS  equations  (balanced  model)  given  in  [51; 

4—  experimental  result  [7]  ; 

5—  adiabatic  wall 
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TGURE  7 (a) ,  Distribution  of  surface  Pressure  in  laminar  flow 
ver  three-dimensional  compression  corner,  with  M 


1 — solution  of  the  simpUfiad  NS  equations  given  in  this  paper 

2—  result  for  the  two-dimensional  case; 

3 —  adiabatic  wall 


Figure  7(b) .  Distribution  of  surface  drag  in  laminar  flow  over 
three-dimensional  compression  corner,  with  =  3. 

1 —  solution  of  the  simplified  NS  equations  given  in  this  paper; 

2 —  result  for  the  two-dimensional  case; 

3 —  adiabatic  wall 
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